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Abstract. We deal with entropy solutions to the Cauchy problem 
U , for the isentropic compressible Euler equations in the space-periodic 

case. In more than one space dimension, the methods developed 
. by De Lellis-Szekelyhidi enable us to show failure of uniqueness 

—j on a finite time-interval for entropy solutions starting from any 

^ — ' continuously differentiable initial density and suitably constructed 

bounded initial linear momenta. 

57 
< 

^ ■ 1. Introduction 

In this note, we deal with the Cauchy Problem for the isentropic com- 
pressible Euler equations in the space-periodic setting. Given any con- 
tinuously differentiable initial density, we can construct bounded initial 
linear momenta for which admissible solutions are not unique in more 

l>" ' than one space dimension. 

'xT ■ 

\ We first introduce the isentropic compressible Euler equations of gas 

dynamics in n space dimensions, n > 2 (cf. Section 3.3 of [3]). They 
£NJ \ are obtained as a simplification of the full compressible Euler equations, 

by assuming the entropy to be constant. The state of the gas will be 
described through the state vector 



JH ; V=(p,m) 



whose components are the density p and the linear momentum m. The 
balance laws in force are for mass and linear momentum. The resulting 
system, which consists of n + 1 equations, takes the form: 



dtp + div x m = 

d t m + div x (^)+V x .b(p)] = 
p(-,0) = p° 
m(-,0) = m 



(1.1) 



The pressure p is a function of p determined from the constitutive ther- 
modynamic relations of the gas in question. A common choice is the 
polytropic pressure law 

P(p) = k P J 
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with constants k > and 7 > 1. The set of admissible values is P = 
{p > 0} (cf. [3] and |19|). The system is hyperbolic if 

p'(p) > 0. 

In addition, thermodynamically admissible processes must also satisfy 
an additional constraint coming from the energy inequality 



1 I |2 

1 \m\ 



dt I pe{p) + 2~~p~) +div2 



. . 1 Iml 2 pip) \ 



< (1.2) 



where the internal energy e : R + — > M is given through the law p(r) = 
r 2 e'(r). The physical region for (11. ip is {(p,m)\ \m\ < for some 

constant > 0. For p > 0, u = m/p represents the velocity of the fluid. 

We will consider, from now on, the case of general pressure laws given 
by a function p on [0, 00 [, that we always assume to be continuously 
differentiable on [0, 00 [. The crucial requirement we impose upon p is that 
it has to be strictly increasing on [0,oo[. Such a condition is meaningful 
from a physical viewpoint since it is a consequence of the principles of 
thermodynamics. 

Now, we will disclose the content of this note. Using some techniques 
introduced by De Lellis-Szekelyhidi (cf. [1] and [5]) we can consider any 
continuously differentiable periodic initial density p° and exhibit suitable 
periodic initial momenta mP for which space-periodic weak admissible 
solutions of (jl.ip are not unique on some finite time-interval. 

Theorem 1.1. Let n > 2. Then, for any given function p and any 
given continuously differentiable periodic initial density p°, there exist a 
bounded periodic initial momentum m° and a positive time T for which 
there are infinitely many space-periodic admissible solutions (p, m) of 
(fTTjl on M. n x [0,T[ with p G C^K" x [0,T[). 

Remark 1.2. Indeed, in order to prove Theorem \l.l\ it would be enough 
to assume that the initial density is a Holder continuous periodic func- 
tion: p° G C°' Q (M n ) (cf Proof of Proposition^^. 



Some connected results are obtained in [5] (cf. Theorem 2 therein) as 
a further consequence of their analysis on the incompressible Euler equa- 
tions. Inspired by their approach, we adapt and apply directly to (II. ip 
the method of convex integration combined with Tartar's programme on 
oscillation phenomena in conservation laws (see |21] and |13j). In this 
way, we can show failure of uniqueness of admissible solutions to the 
compressible Euler equations starting from any given continuously dif- 
ferentiable initial density. For a survey on these /i-principle-type results 
in fluid dynamics we refer the reader to [6]. 
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The paper is organised as follows. Section 2 is an overview on the 
definitions of weak and admissible solutions to (jl.ip and gives a first 
glimpse on how Theorem 11.11 is achieved. Section 3 is devoted to the 
reformulation of a simplified version of the isentropic compressible Euler 
equations as a differential inclusion and to the corresponding geometrical 
analysis. In Section 4 we state and prove a criterion (Proposition [JTJ to 
select initial momenta allowing for infinitely many solutions. The proof 
builds upon a refined version of the Baire category method for differential 
inclusions developed in [5] and aimed at yielding weakly continuous in 
time solutions. Section 5 and 6 contain the proofs of the main tools 
used to prove Proposition 14.11 In Section 7, we show initial momenta 
satisfying the requirements of Proposition 14.11 Finally, in Section 8 we 
prove Theorem II .11 by applying Proposition 14.11 

2. Weak and admissible solutions to 
the isentropic euler system 

The deceivingly simple-looking system of first-order partial differential 
equations (jl.ip has a long history of important contributions over more 
than two centuries. We recall a few classical facts on this system (see for 
instance [3] for more details). 

• If p° and m° are "smooth" enough (see Theorem 5.3.1 in [3]), 
there exists a maximal time interval [0, T[ on which there exists 
a unique "smooth" solution (p,m) of (jl.ip (for < t < T). In 
addition, if T < oo, and this is the case in general, (p, m) becomes 
discontinuous as t goes to T. 

• If we allow for discontinuous solutions, i.e., for instance, solu- 
tions (p,m) E L°° satisfying (jl.ip in the sense of distributions, 
then solutions are neither unique nor stable. More precisely, one 
can exhibit sequences of such solutions which converge weakly in 
L°° — * to functions which do not satisfy (jl.ip . 

• In order to restore the stability of solutions and (possibly) the 
uniqueness, one may and should impose further restrictions on 
bounded solutions of (jl.ip . restrictions which are known as (Lax) 
entropy inequalities. 

This note stems from the problem of better understanding the effi- 
ciency of entropy inequalities as selection criteria among weak solutions. 

Here, we have chosen to emphasize the case of the flow with space 
periodic boundary conditions. For space periodic flows we assume that 
the fluid fills the entire space lR n but with the condition that m, p are 
periodic functions of the space variable. The space periodic case is not 
a physically achievable one, but it is relevant on the physical side as a 
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model for some flows. On the mathematical side, it retains the com- 
plexities due to the nonlinear terms (introduced by the kinematics) and 
therefore it includes many of the difficulties encountered in the general 
case. However the former is simpler to treat because of the absence of 
boundaries. Furthermore, using Fourier transform as a tool simplifies 
the analysis. 

Let Q = [0, If, n > 2 be the unit cube in R n . We denote by H™(Q), 
m E N, the space of functions which are in H^iW 1 ) and which are 
periodic with period Q: 

m{x + /) = m(x) for a.e. x E W 1 and every I E Z n . 

For m = 0, H®(Q) coincides simply with L 2 (Q). Analogously, for every 
functional space X we define X p (Q) to be the space of functions which 
are locally (over M n ) in X and are periodic of period Q. The functions 
in H p n (Q) are easily characterized by their Fourier series expansion 

H™(Q) = lme L 2 p {Q) : ^ \k\ 2m \m(k)\ 2 < oo and m(0) = 

(2-1) 

where m : Z n — > C n denotes the Fourier transform of m. We will use 
the notation H(Q) for H p {Q) and H W (Q) for the space H(Q) endowed 
with the weak L 2 topology. 

Let T be a fixed positive time. By a weak solution of (jl.ip on lR n x [0, T[ 
we mean a pair (p,m) E L°°([0,T[; L2°(Q)) satisfying 

|m(x,i)| < Rp(x,t) for a.e. (x,f) £l"x [0,T[ and some i? > 0, 

(2.2) 

and such that the following identities hold for every test functions if) E 
C-([0,T[;C-(g)), 0ECr([O,T[;C-(Q)): 

f [ [pd t i> + m-X7 x ip}dxdt+ [ p°(x)ip(x,0)dx = (2.3) 
^0 JQ JQ 



'0 JQ 



TYli Q^) TXXi 

( — - — , V x <t> ) + p(p) div a 



dxdt 



+ / m°(x) ■ <j)(x,Q)dx = 0. 



(2.4) 



For n > 2 the only non-trivial entropy is the total energy n = pe(p) + 
h^~p~ wn i cn corresponds to the flux ^ = \s(p) + — h 



m. 
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Then a bounded weak solution (p, m) of (jl.ip satisfying (jl.2p in the 
sense of distributions, i.e. satisfying the following inequality 



T 



pe(p) + I 9t<p + (e(p) + iM- + ^ J m • V x y 



Q 




pW) + o J — L ^ °- ( 2 - 5 ) 



for every nonnegative 99 € C£°([0, T[; Cf°(Q)), is said to be an entropy 
(or admissible) solution of f 1 1 . 1 j) . 

The lack of entropies is one of the essential reasons for a very lim- 
ited understanding of compressible Euler equations in dimensions greater 
than or equal to 2. 

A recent paper by De Lellis-Szekelyhidi gives an example in favour 
of the conjecture that entropy solutions to the multi-dimensional com- 
pressible Euler equations are in general not unique. Showing that this 
conjecture is true has far-reaching consequences. The entropy condition 
is not sufficient as a selection principle for physical/unique solutions. 
The non-uniqueness result by De Lellis-Szekelyhidi is a byproduct of 
their new analysis of the incompressible Euler equations based on its 
formulation as a differential inclusion. They first show that, for some 
bounded compactly supported initial data, none of the classical admis- 
sibility criteria singles out a unique solution to the Cauchy problem for 
the incompressible Euler equations. As a consequence, by constructing 
a piecewise constant in space and independent of time density p, they 
look at the compressible isentropic system as a "piecewise incompress- 
ible" system (i.e. still incompressible in the support of the velocity field) 
and thereby exploit the result for the incompressible Euler equations to 
exhibit bounded initial density and bounded compactly supported initial 
momenta for which admissible solutions of (jl.ip are not unique (in more 
than one space dimension). 

Inspired by their techniques, we give a further counterexample to the 
well-posedeness of entropy solutions to (II. ip . Our result differs in two 
main aspects: here the initial density can be any given "regular" function 
and remains "regular" forward in time while in [5] the density allowing 
for infinitely many admissible solutions must be chosen as piecewise con- 
stant in space; on the other hand we are not able to deal with compactly 
supported momenta (indeed we work in the periodic setting), hence our 
non-unique entropy solutions are only locally I? in contrast with the 
global-L 2 -in-space property of solutions obtained in [5]. Moreover, we 
have chosen to study the case of the flow in a cube of M. n with space 
periodic boundary conditions. This case leads to many technical simpli- 
fications while retaining the main structure of the problem. 
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More precisely, we are able to analyze the compressible Euler equa- 
tions in the framework of convex integration. This method works well 
with systems of nonlinear PDEs such that the convex envelope (in an 
appropriate sense) of each small domain of the submanifold represent- 
ing the PDE in the jet-space (see [8] for more details) is big enough. 
In our case, we consider a simplification of system (jl.ip . namely the 
semi-stationary associated problem, whose submanifold allows a convex 
integration approach leading us to recover the result of Theorem 11.11 

We are interested in the semi-stationary Cauchy problem associated 
with the isentropic Euler equations (simply set to the time derivative 
of the density in (jl.ip and drop the initial condition for p) : 

div x m = 

d t m + div x (ass) + VMP)} = (2.6) 
m(-,0) = m . 

A pair (p,m) G L™(Q) x L°°([0, T[; L™(Q)) is a weak solution on R n x 
[0, T[ of (|2.6|) if m(-, t) is weakly-divergence free for almost every < t < 
T and satisfies the following bound 

\m(x,t)\ < Rp(x) for a.e. (x,t) £l"x [0, T[ and some R > 0, 

(2.7) 

and if the following identity holds for every <f) £ C£°([0,T[; C™(Q)): 

cT 



m-d t <p+ { — - — , V x (f> ) + p(p) div a 



dxdt 



o JQ 

+ / m°(x) ■ <f){x,0)dx = 0. (2.8) 
JQ 

A general observation suggests us that a non-uniqueness result for 
weak solutions of (|2.6p whose momentum's magnitude satisfies some 
suitable constraint could lead us to a non-uniqueness result for entropy 
solutions of the isentropic Euler equations (jl.ip . Indeed, the entropy 
solutions we construct in Theorem 11.11 come from some weak solutions 
of (1^61). 



Theorem 2.1. Let n > 2. Then, for any given function p, any given 
density po E Cp(Q) and any given finite positive time T, there exists a 
bounded initial momentum m° for which there are infinitely many weak 
solutions (p,m) G C£(Q) x C([0, T]; H W {Q)) of onW 1 x [0,T[ with 
density p(x) = po(x). 

In particular, the obtained weak solutions m satisfy 

\m(x,t)\ 2 = p {x)x(t) a.e. inR n x [0,T[, (2.9) 
\m°(x)\ 2 = p (x)x(0) a.e. inR", (2.10) 

for some smooth function \- 
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An easy computation shows how, by properly choosing the function 
X in (l2.9D - (l2.10p . the solutions (po,m) of fj2.6[) obtained in Theorem 12.11 
satisfy the admissibility condition ()2.5[) . 

Theorem 2.2. Under the same assumptions of Theorem \2.1\ there ex- 
ists a maximal time T > such that the weak solutions (p, m) of (|2.6I) 
( coming from Theorem \2. 1\) satisfy the admissibility condition (12. 5ft on 



Our construction yields initial data m° for which the nonuniqueness 
result of Theorem 1 1 . 1 1 holds on any time interval [0, T[, with T < T. How- 
ever, as pointed out before, for sufficiently regular initial data, classical 
results give the local uniqueness of smooth solutions. Thus, a fortiori, 
the initial momenta considered in our examples have necessarily a certain 
degree of irregularity. 



This section is devoted to a qualitative analysis of the isentropic com- 
pressible Euler equations in a semi-stationary regime (i.e. (|2.6p ). 

As in [3] we will interpret the system (|2.6|) in terms of a differential 
inclusion, so that it can be studied in the framework combining the plane 
wave analysis of Tartar, the convex integration of Gromov and the Baire's 
arguments. 

3.1. Differential inclusion. The system (12.6j) can indeed be naturally 
expressed as a linear system of partial differential equations coupled with 
a pointwise nonlinear constraint, usually called differential inclusion. 

The following Lemma, based on Lemma 2 in [5], gives such a reformu- 
lation. We will denote by S n the space of symmetric n x n matrices, by 
Sq the subspace of S n of matrices with null trace, and by I n the n x n 
identity matrix. 



Lemma 3.1. Letm G L°°([0, T\\ L~(Q; R n )), U G L°°([0, T]; L|°(Q; «S£)) 
and q G L°°([0,T]; L°°(Q;R+)) such that 



If (m, U,q) solve (|3.ip and in addition there exists p G L^°(R n ;]R + ) such 
that (|2.7p holds and 



[0,T[. 



3. Geometrical analysis 



diyj; m = 



dtin + diva U + V x q = 0. 



(3.1) 



2 



U = 



m 



I n a.e. in R n x [0,T] 



P 



np 



2 



q = p(p) + 



m 



a.e. in R n x [0,T] 



(3.2) 



np 
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then m and p solve (|2,6p distributionally. Conversely, if m and p are 

T- 1 € I nand q = P (p) + ^ 



I |2 II 2 

weak solutions of (|2,6p . then m, U = m ® m — ^-I n and q = p(p) + 



solve JO-G321. 




In Lemma 13.11 we made clear the distinction between the augmented 
system (|3.ip . whose linearity allows a plane wave analysis, and the non- 
linear pointwise constraint (I3.2p . which leads us to study the graph below. 

For any given p e]0,oo[, we define the following graph 

5 "xI+:[J= - ^Ll n , 

p np 

(3.3) 

The key of the forthcoming analysis is the behaviour of the graph K p 
with respect to the wave vectors associated with the linear system (|3.ip : 
are differential and algebraic constraints in some sense compatible? 

For our purposes, it is convenient to consider "slices" of the graph K p , 
by considering vectors m whose modulus is subject to some /^-depending 
condition. Thus, for any given \ £ ]R + , we define: 



K px := \(m,U,q) 6 l n x iS n x ]R + : U = - ^I n , 

p np 

I |2 1 

<7=p(p)H , m =px>- (3.4) 

np 



3.2. Wave cone. Following Tartar's framework |21j . we consider a sys- 
tem of first order linear PDEs 

Y,AidiZ = (3.5) 

i 

where z is a vector valued function and the A, are matrices. Then, 
planewave solutions to (13. 5p are solutions of the form 

z(x) = a/i(x • ^), (3.6) 

with /i : M — > R. In order to find such solutions, we have to solve 
the relation Y2i £iAia = 0, where £j is the oscillation frequency in the 
direction i. The set of directions a for which a solution £ 7^ exists 
is called wane cone A of the system (|3.5p : equivalently A characterizes 
the directions of one dimensional high frequency oscillations compatible 
with ([33]) . 
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The system (|3.ip can be analyzed in this framework. Consider the 
(n + 1) x (n + 1) symmetric matrix in block form 



Note that, with the new coordinates y = (x,t) G M n+1 , the system (|3.ip 
can be easily rewritten as div y M = (the divergence of M in space-time 
is zero). Thus, the wave cone associated with the system (|3.ip is equal 
to 



Indeed, the relation Y2i = f° r the system (|3.ip reads simply 

as M ■ (£, c) = 0, where (£, c) G R n x M (£ is the space- frequency and 
c the time- frequency) : this equation admits a non-trivial solution if M 
has null determinant, hence (|3.8p . 

3.3. Convex hull and geometric setup. Given a cone A, we say that 
K is convex with respect to A if, for any two points A,B G K with 
B — A G A, the whole segment [A, B] belongs to K. The A-convex 
hull of K P:X is the smallest A-convex set containing K PtX , i.e. the 

set of states obtained by mixture of states of K PiX through oscillations 
in A-directions (Gromov who works in the more general setting of 
jet bundles, calls this the P- convex hull). The key point in Gromov's 
method of convex integration (which is a far reaching generalization of 
the work of Nash |17| and Kuiper |14j on isometric immersions) is that 
(|3.5p coupled with a pointwise nonlinear constraint of the form z G K 
a.e. admits many interesting solutions provided that the A-convex hull 
of K, i^ A , is sufficiently large. In applications to elliptic and parabolic 
systems we always have K A = K so that Gromov's approach does not 
directly apply. For other applications to partial differential equations it 
turns out that one can work with the A-convex hull defined by duality. 
More precisely, a point does not belong to the A- convex hull defined by 
duality if and only if there exists a A-convex function which separates it 
from K. A crucial fact is that the second notion is much weaker. This 
surprising fact is illustrated in |13| . 

In our case, the wave cone is quite large, therefore it is sufficient to 
consider the stronger notion of A-convex hull, indeed it coincides with 
the whole convex hull of K PtX . 



Lemma 3.2. For any S G S n let X ma x(S) denote the largest eigenvalue 
of S. For (p, m, U) G M+ x R n x S™ let 




(3.7) 




(3.8) 




(3.9) 
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Then, for any given p, x G R + , the following holds 



(i) e(p, •, •) : R n x Sq — > R is convex; 

(ii) J=£ < e(p, m, 17), eguaZtfj, i/ and on/y £/*7 = » - ^7 n ; 

(iii) |£7| < (n — l)e(p,m,U), with \U\ being the operator norm of 
the matrix; 

(iv) the ^-sublevel set of e defines the convex hull of K PtX , i.e. 

K c ° x = | (m, U, q) G R n x <S£ x R + : e(p, m,U)<^, 

q = p(p) + -\ (3.10) 
n I 

andK PiX = K™ x n{\m\ 2 = PX }. 

For the proof of (i)-(iv) we refer the reader to the proof of Lemma 3.2 
in [5]: the arguments there can be easily adapted to our case. 

We observe that, for any £ R + , the convex hull Kp° x lives in 
the hyperplane H of W 1 x S™ x R+ defined by H := |(m, U, q) G M n x 

Sq x R + : q = p{p) + %\. Therefore, the interior of Kp° x as a subset 
of R n x Sq x R + is empty. This seems to prevent us from working in 
the classical framework of convex integration, but we can overcome this 
apparent obstacle. 

For any p, \ G R + , we define the hyperinterior of Kp° x , and we denote 
it with "hint Kp° x \ as the following set 

hint K c ° := J (m, U, q) G R n x x R+ : e(p, m,U) <-, 

n 

q = p(p) + -\. (3.11) 
n 



In the framework of convex integration, the larger the A-convex hull 
of K P)X is, the bigger the breathing space will be. How to "quantify" the 
meaning of a "large" A-convex hull in our context? The previous defini- 
tion provides an answer: the A-convex hull of K PtX will be "large" if its 
hyperinterior is nonempty. The wave cone of the semi-stationary Euler 
isentropic system is wide enough to ensure that the A-convex hull of K p<x 
coincides with the convex hull of K„ x and has a nonempty hyperinterior. 
As a consequence, we can construct irregular solutions oscillating along 
any fixed direction. For our purposes, it will be convenient to restrict to 
some special directions in A, consisting of matrices of rank 2, which are 
not stationary in time, but are associated with a constant pressure. 
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Lemma 3.3. Let c,d£ W 1 ' with \c\ = \d\ and c / d, and let p G M + . 
Then (c-d,^-^,0) G A. 

Proof. Since the vector + d, — ^ l c l + c i s j n the kernel of the ma- 
trix 

P P 

c-d J 

C has indeed determinant zero, hence — d, ^ — 0^ G A. □ 

Now, we introduce some important tools: they allow us to prove that 
Kp X = Kp° x ls sufficiently large, thus providing us room to find many 
solutions for (|3~TT)-(|3~21). 



As first, we define the admissible segments as segments in W 1 x5q xR + 
whose directions belong to the wave cone A for the linear system of PDEs 
(|3.ip and are indeed special directions in the sense specified by Lemma 
1331 



Definition 3.4. Given p, \ G ^ + we °~ an admissible segment for 
(Pj x) if c is a line segment in W l x 5J x M + satisfying the following 
conditions: 

• a is contained in the hyperinterior of K c p ° x ; 

• a is parallel to (c- d,^ - ^,6) for some c,d G R n with 
\c\ 2 = \d\ 2 = px and c ^ ±d. 

The admissible segments defined above correspond to suitable plane- 
wave solutions of (|3.ip . The following Lemma ensures that, for any 
p, X G M + , the hyperinterior of Kp° x is " sufficiently round " with respect 
to the special directions: given any point in the hyperinterior of Kp° x , 
it can be seen as the midpoint of a sufficiently large admissible segment 
for (p,x). 

Lemma 3.5. There exists a constant F = F(n) > such that for any 
p, X G R + and for any z = (m, U, q) G hint Kp° x there exists an admissi- 
ble line segment for (p, x) 

a= [(m,U,q) - (m,U,0),{m,U,q) + (m,U,0)) (3.12) 

such that 

F 

The proof rests on a clever application of Caratheodory's theorem for 
convex sets and can be carried out, with minor modifications, as in [5] 
(cf. Lemma 6 therein). 

As an easy consequence of the previous Lemma, we can finally estab- 

f PjX coincides with Kp° x . 



F ( sl\ 

m\ > I px — \m\ 
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Proposition 3.6. For all given p, x G the A-convex hull of K PtX 
coincides with the convex hull of K PiX . 

Proof. Recall that, given p, x £ K + , we denote the A-convex hull of K p ^ x 
with K p x . Of course K px C K p ° x , hence we have to prove the opposite 
inclusion, i.e. K p ° x C K px . For every z € K p ° x we can follow the 
procedure in the proof of Lemma [3.5l (cf. [5]) and write it as z = ^ • XjZj, 
with (zj)i<j<N+i in A AX , (Xj)i<j<N+i in [0,1] and Aj = 1- Again, 
we can assume that Ai = max-,- A-,-. In case Ai = 1 then z = Z\ € i^p,x C 
and we can already conclude. Otherwise (i.e. when Ai G (0, 1)) we 
can argue as in Lemma [3 .51 so to find an admissible segment a for (p, %) of 
the form (j3. 12[) . Since we aim at writing z as a A-barycenter of elements 
of K px , we "play" with these admissible segments by prolongations and 
iterative constructions until we get segments with extremes lying in K PtX . 
More precisely: we extend the segment a until we meet d hint-R^" thus 
obtaining z as the barycenter of two points (wo, w±) with (wq — W\) £ A 
and such that every ujj = (m^, C/j, cjj), i = 0, 1, satisfies either |m.j| 2 = p% 
or \m.i\ 2 < px and e(p, nii, Ui) = x/ n - 

In the first case, £/, - ( m '® m » - ^^I n ) > 0, and since it is a null- 
trace-matrix it is identically zero, whence W{ E Kp,x ( n °te that in the 
construction of a the g-direction remains unchanged, hence q± = p(p) + 

i)- 

In the second case, i.e. when \mi\ < px and e(p,mi,Ui) = xl n i 
we apply again Lemma 13.51 and a limit procedure to express Wi as 
barycentre of (to* Oj with (io^o — ^i,l) £ A and such that every 
UJi.fc = {rm^,Ui^,qi,k), k = 0,1, will satisfy either |mj jfc | 2 = or 
A2(p, rrii,k, U it k) = e(p, m i)fc , ET^) = x/n, where Ai(p, m, E7) > A 2 (p, m, U) > 

> X n (p,m,U) denote the ordered eigenvalues of the matrix m ® m —JJ 

(note that Xi(p,m,U) = e(p, m, [/)). Now, we iterate this procedure of 
constructing suitable admissible segments for (p, %) until we have writ- 
ten z as A-barycenter of points (m,U,q) satisfying either |m| 2 = px or 
X n (p,m, U) = x/ n an d therefore all belonging to K P)X as desired. □ 



4. A CRITERION FOR THE EXISTENCE OF INFINITELY MANY 

SOLUTIONS 

The following Proposition provides a criterion to recognize initial data 
m° which allow for many weak admissible solutions to f)l . 1 [) . Its proof 
relies deeply on the geometrical analysis carried out in Section 3. The 
underlying idea comes from convex integration. The general principle of 
this method, developed for partial differential equations by Gromov [11] 
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and for ordinary differential equations by Filippov |10] . consists in the 
following steps: given a nonlinear equation £(z), 

• (i) we rewrite it as (C(z) A z € K) where £ is a linear equation; 

• (ii) we introduce a strict subsolution zq of the system, i.e. satis- 
fying a relaxed system (C(zq) A z S U); 

• (hi) we construct a sequence (zk)k<=N approaching K but staying 
in U\ 

• (iv) we pass to the limit, possibly modifying the sequence (zk) in 
order to ensure a suitable convergence. 

Step (i) has already been done in Section 3.1. The choice of zq will be 
specified in Sections 7-8. Here, we define the notion of subsolution for 
an appropriate set U, we construct an improving sequence and we pass 
to the limit. The way how we construct the approximating sequence will 
be described in Section 6 using some tools from Section 5. 

One crucial step in convex integration is the passage from open sets 
K to general sets. This can be done in different ways, e.g. by the Baire 
category theorem (cf. |18j). a refinement of it using Baire- 1 functions or 
the Banach-Mazur game |12| or by direct construction |20] , Whatever 
approach one uses the basic theme is the same: at each step of the 
construction one adds a highly oscillatory correction whose frequency 
is much larger and whose amplitude is much smaller than those of the 
previous corrections. 

In this section, we achieve our goals following some Baire category 
arguments as in pQ: they are morally close to the methods developed by 
Bressan and Flores in [T) and by Kirchheim in |12| . 

In our framework the initial data will be constructed starting from 
solutions to the convexified (or relaxed) problem associated to (|2,6p . i.e. 
solutions to the linearized system (|3.ip satisfying a "relaxed" nonlinear 
constraint (I3.2P (i.e. belonging to the hyperinterior of the convex hull of 
the "constraint set"), which we will call subsolutions. 

As in [3], our application shows that the Baire theory is comparable in 
terms of results to the method of convex integration and they have many 
similarities: they are both based on an approximation approach to tackle 
problems while the difference lies only in the limit arguments, i.e. on the 
way the exact solution is obtained from better and better approximate 
ones. These similarities are clarified by Kirchheim in |12| . where the 
continuity points of a first category Baire function are considered; a 
comparison between the two methods is drawn by Sychev in [20J. 

Here, the topological reasoning of Baire theory is preferred to the it- 
eration technique of convex integration, since the first has the advantage 
to provide us directly with infinitely many different solutions. 
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Proposition 4.1. Let po G Cp(Q; M + ) be a given density function and let 
T be any finite positive time. Assume there exist (mo,Uo,qo) continuous 
space-periodic solutions of (|3.ip on IR n x]0, T[ with 

m G C([0,T];H w (Q)), (4.1) 
and a function % G C°°([0, T]; K + ) such that 

e(p (x),m (x,t),U (x,t)) < ^ /or a// (rr, t) G M n x]0, T[, (4.2) 

n 

qo(x,t)= P (po(x)) + ^ for all (x,t) eR n x]0,T[. (4.3) 
n 



Then there exist infinitely many weak solutions (p, m) of the system (|2.6p 
in W 1 x [0, T[ with density p{x) = po(x) and such that 

meC([0,T];H w (Q)), (4.4) 
m(-,t) = m (-,t) for t = 0,T and for a.e. x el", (4.5) 
\m(x,t)\ 2 = p (x)x{t) for a.e. (x,t) G M n x]0,T[. (4.6) 

4.1. The space of subsolutions. We define the space of subsolutions 
as follows. Let po and x be given as in the assumptions of Proposi- 
tion 14.11 Let mo be a vector field as in Proposition 14.11 with associ- 
ated modified pressure qo and consider space-periodic momentum fields 
m :M. n x [0, T] —?■ I" which satisfy 

divm = 0, (4.7) 

the initial and boundary conditions 

m(x,0) =m (x,0), (4.8) 
m(x,T) = m (x,T), (4.9) 

(4.10) 

and such that there exists a continuous space-periodic matrix field U : 
R n x]0,T[-^ S% with 

e(p (x),m(x,t),U(x,t)) < - for all (x,t) G R n x]0,T[, 

n (4.11) 
d t m + div x U + V x q = in R n x [0, T]. 

Definition 4.2. Let Xq be the set of such linear momentum fields, i.e. 
X = |m G C°(]0, T[; Cj(Q)) fl C([0, T]; fl^Q)) : 

(|4TT|) - ([4~TLj) are safes/iedj (4.12) 

and let X be the closure of Xq in C([0, T]; H W (Q). Then Xq will be the 
space of strict subsolutions. 
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As po is continuous and periodic on R n and \ is smooth on [0, 7], 
there exists a constant G such that x(t) fn Po(x)dx < G for all t G [0, 7]. 
Since for any m G Xq with associated matrix field U we have that (see 
Lemma [3.21 - (ii)) 



\m(x,t)\ dx < npo(x)e(po(x),m(x,t),U(x,t))dx 
Q JQ 

< X{t) / po{x)dx for all t G [0,7], 
JQ 

we can observe that Xq consists of functions m : [0, T] — > H(Q) taking 
values in a bounded subset B of H(Q). Without loss of generality, we 
can assume that B is weakly closed. Then, B in its weak topology is 
metrizable and, if we let dg be a metric on B inducing the weak topology, 
we have that (B,ds) is a compact metric space. Moreover, we can define 
on Y := C([0, T], (B,dB)) a metric d naturally induced by ds via 

d(h,f 2 ):= maxd B (fi(-,t),f 2 (.-,t)). (4.13) 
*e[o,T] 

Note that the topology induced on Y by d is equivalent to the topology of 
Y as a subset of C([0, 7]; H w ). In addition, the space (Y, d) is complete. 
Finally, X is the closure in (Y, d) of Xq and hence (X, d) is as well a 
complete metric space. 

Lemma 4.3. If m € X is such that \m(x,t)\ 2 = pQ(x)x(t) for almost 



every (x,t) G M n x]0,T[, then the pair (pQ,m) is a weak solution of fj2 .6|) 
in R n x [0,T[ satisfying fO]> - (|Q]l - (|4T5|> . 

Proof. Let m£lbe such that |m(j;,t)| 2 = po^X^) f° r almost every 
(x,t) € M n x]0,T[. By density of Xq, there exists a sequence {m^} C ATo 
such that mfc -4 m in X. For any G Ao let Uk be the associated 
smooth matrix field enjoying (14. lip . Thanks to Lemma 13.21 (iii) and 
(|4,lip . the following pointwise estimate holds for the sequence {Uk} 

(ft — \\ — v 

|£4L < {n-l)e{pQ,m k ,U k ) < ± 

As a consequence, {Uk} is uniformly bounded in L°°([0, 7]; L^°(Q)); by 
possibly extracting a subsequence, we have that 

Uk^U in L°°([0,T];£~(Q)). 



Note that hintETp° = Kp° tX is a convex and compact set by Lemma 
3.2l -(i)-(ii)-(iii). Hence, m £ X with associated matrix field U solves 
(13. ip on R n x [0, 7] for q = qo and (m, ?7, go) takes values in Kp° x almost 
everywhere. If, in addition, \m(x,t)\ 2 = pQ(x)x(t), then (rn, U, qo)(x, t) G 
ifp jX a.e. in M n x [0,7] (cf. Lemma l3.2I (iv)). Lemma 13.11 allows us to 
conclude that (pQ,m) is a weak solution of (|2,6p in M n x [0, T[. Finally, 
since mj; — > m in C([0, T]; H W (Q)) and |m(x,t)| 2 = po(^)x(0 f° r almost 
every (x, i) G M n x]0, T[, we see that m satisfies also (|4T4) - (|43]) - (|Q]) . □ 
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Now, we will argue as in [1] exploiting Baire category techniques to 
combine weak and strong convergence (see also |12]). 

Lemma 4.4. The identity map I : (X,d) — > L 2 ([0,T]; H(Q)) defined by 
m —7- m is a Baire-1 map, and therefore the set of points of continuity is 
residual in (X, d). 

Proof. Let <f) r (x,t) = r~( n+1 ^(/)(rx,rt) be any regular spacetime convo- 
lution kernel. For each fixed m G X, we have 

<fi r * m — > m strongly in L 2 (H) as r — > 0. 

On the other hand, for each r > and mk G X , 

nrik —> rn implies (j) r * m/% — > (j) r * m in L 2 (H). 

Therefore, each map I r : (X, d) — > L 2 (H), m — > (j) r *m is continuous, and 
I(m) = lim r _>o I r {m) for all m £ X. This shows that / : (X, d) — > L 2 {H) 
is a pointwise limit of continuous maps; hence it is a Baire-1 map. As a 
consequence, the set of points of continuity of / is residual in (X, d) (cf. 

M)- □ 



4.2. Proof of Proposition 14.11 We aim to show that all points of 
continuity of the identity map correspond to solutions of (|2.6p enjoying 
the requirements of Proposition 14.11 Lemma 14.41 will then allow us to 
prove Proposition 14. II once we know that the cardinality of X is infinite. 
In light of Lemma 14.31 for our purposes it suffices to prove the following 
claim: 



CLAIM. If m € X is a point of continuity of /, then 

\m(x,t)\ 2 = po(x)x(t) for almost every (x,t) E R"x]0,T[. (4.14) 

□ 



Note that proving (|4.14p is equivalent to prove that || w|| i 2(g x r T n = 

( T \ 1 / 2 2 

( Iq Jo Po{ x )x(t)dtdx) , since for any m G X we have |m(a;,t)| < 

Po(x)x(t) for almost all (x,t) G K™ x [0,T]. Thanks to this remark, the 

claim is reduced to the following lemma (cf. Lemma 4.6 in |3]), which 

provides a strategy to move towards the boundary of Xq. given m G Xq, 

we will be able to approach it with a sequence inside Xq but closer than 

m to the boundary of Xq. 



Lemma 4.5. Let po,x be given functions as in Proposition ^. 1\ Then, 
there exists a constant (3 = (3(n) such that, given m G Xq, there exists a 
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sequence {m^} C Xq with the following properties 

ll ll 2 >, ll ll 2 

H m fcHL 2 (Qx[0,T]) - W m \\L 2 (Qx[0,T}) 

+ P (J po(x)x(t)dtdx - ||m||| 2(Qx[0iT]) 

(4.15) 

and 

m k ^ m inC([0,T],H w (Q)). (4.16) 



The proof is postponed to Section [6j Let us show how Lemma 14.51 
implies the claim. As in the claim, assume that m G X is a point of 
continuity of the identity map /. Let {m^} C Xq be a fixed sequence that 
converges to m in C([0, T], H W (Q)). Using Lemma 14.51 and a standard 
diagonal argument, we can find a second sequence {fhk} yet converging 
to m in X and satisfying 

^minf \\™>k\\l*(Qx[o,T]) > !™in f ^ H m fcllL2(Qx[o,r]) 

+ P( [ [ Po( x )x(t)dtdx - ||m fe ||^2 ( Q x[0)T] ) 

According to the hypothesis, / is continuous at m, therefore both m*, 
and rhk converge strongly to m and 

2 



\ m \\L2(Qx[0,T}) ^ ll m llL 2 (Qx[0,T]) 



+ /3 (/ lo P0 ^ X ^ dtdX ~ H m lli 2 (Qx[0,T])) 



2 



Hence ||m-||i2(Q X [ t]) = {^Iq Iq Po( x )x{t)dtdx^ and the claim holds 
true. Finally, since the assumptions of Proposition 14. 1 1 ensure that Xq is 
nonempty, by Lemma [4.51 we can see that the cardinality of X is infinite 
whence the cardinality of any residual set in X is infinite. In particular, 
the set of continuity points of / is infinite: this and the claim conclude 
the proof of Proposition 14.11 



1/2 



5. Localized oscillating solutions 



The wild solutions are made by adding one dimensional oscillating 
functions in different directions A £ A. For that it is needed to localize 
the waves. More precisely, the proof of Lemma 14.51 relies on the con- 
struction of solutions to the linear system (|3.ip . localized in space-time 
and oscillating between two states in along a given special direc- 

tion A € A. Aiming at compactly supported solutions, one faces the 
problem of localizing vector valued functions: this is bypassed thanks 
to the construction of a "localizing" potential for the conservation laws 
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(|3.ip . This approach is inherited from [5]. As in [3] it could be realized 
for every A € A, but in our framework it is convenient to restrict only to 
special A-directions (cf. [5]): this restriction will allow us to localize the 
oscillations at constant pressure. 

Why oscillations at constant pressure are meaningful for us and needed 
in the proof of Lemma 14.51 / 

Owing to Section 3, in the variables y = (x,t) G R n+1 , the system 
(|3.ip is equivalent to div y M = 0, where M G S n+1 is defined via the 
linear map 



More precisely, this map builds an identification between the set of solu- 
tions (m, U, q) to (|3.ip and the set of symmetric (n+ 1) x (n+ 1) matrices 
M with M (n+1 )( n+1 ) = and tr(M) = q. 

Therefore, solutions of (|3.ip with q = correspond to matrix fields 
M : R n+1 -> ]R(«+i)x(«+i) such that 

di V?/ M = 0, M T = M, M (n+1)(n+1) = 0, tr(M) = 0. (5.2) 

Moreover, given a density p and two states (c,U c ,q c ), {d,Ud,qd) G -Kp 
with non collinear momentum vector fields c and <i having same magni- 
tude (|c| = |(i|), and hence same pressure (q c = q^), then the correspond- 
ing matrices M c and have the following form 



Finally note that tr(M c — M^) = and M c — G A corresponds to a 
special direction. 

The following Proposition provides a potential for solutions of ( 13. ID 
oscillating between two states M c and at constant pressure. It is an 
easy adaptation to our framework of Proposition 4 in [5|. 

Proposition 5.1. Let c,d£ R n such that \c\ = \d\ and c ^ d. Let also 
p £ R. Then there exists a matrix- valued, constant coefficient, homoge- 
neous linear differential operator of order 3 



A{3) : C c °°(R n+1 ) -> C7 c °°(M n+1 ;R( n+1 ) x ( n+1 )) 

such that M = A(d)4> satisfies fl52J) for all (f> G C£°(M n+1 ). Moreover 
there exists rj G R n+1 suc/i £/iai 

• r] is not parallel to e n+ \; 




(5.1) 




and satisfy 
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• if <t>{y) = i>(y • v)> then 

A{d)4>{y) = {M c -M d )4?"{yn). 

We also report Lemma 7 from [5]: it ensures that the oscillations of 
the planewaves generated in proposition 15. II have a certain size in terms 
of an appropriate norm-type-functional. 

Lemma 5.2. Let n 6 M n+1 be a vector which is not parallel to e n+ \. 
Then for any bounded open set B C W 1 

lim / sin 2 (A r r/ • (x, t))dx = — \B\ 

uniformly in t £ R. 

For the proof we refer the reader to [5|. 

6. The improvement step 

We are now about to prove one of the cornerstones of the costruc- 
tion. Before moving forward, let us resume the plan. We have already 
identified a relaxed problem by introducing subsolutions. Then, we have 
proved a sort of " /i-principle" (even if there is no homotopy here) accord- 
ing to which, the space of subsolutions can be "reduced" to the space 
of solutions or, equivalently, the typical (in Baire's sense) subsolution 
is a solution. Once assumed that a subsolution exists, the proof of our 
"/i-principle" builds upon Lemma 14.51 combined with Baire category ar- 
guments. Indeed, we could also prove Proposition 14 . 1 1 by applying itera- 
tively Lemma 14.51 and thus constructing a converging sequence of subso- 
lutions approaching K PtX : this would correspond to the convex integra- 
tion approach. So two steps are left in order to conclude our argument: 
showing the existence of a "starting" subsolution and prove Lemma [4.51 

This section is devoted to the second task, the proof of Lemma 14.51 
while in next section we will exhibit a "concrete" subsolution. 

What follows will be quite technical, therefore we first would like to 
recall the plan: we will add fast oscillations in allowed directions so to 
let |m| 2 increase in average. The proof is inspired by |l]-[5]- 

Proof. [Proof of Lemma I4.5| Let us fix the domain f2 := Q x [0,T]. 
We look for a sequence {m^} C Xq, with associated matrix fields {U^}, 
which improves m in the sense of ([4.15P and has the form 

(m k , U k ) = (m, U) + ^(mjtj, U k ,j) (6.1) 

3 

where every z k j = (fh k j,U k j) is compactly supported in some suitable 
ball B k j(x k j,t k j) C We proceed as follows. 
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Step 1. Let m G Xq with associated matrix field U. By Lemma \3. 51 for 
any(x,i) G £1 we can find a line segment cr^^) := [(m(x,t),U(x,t),qo(x))- 
\x,t), {m(x,t),U{x,t),q (x))+\( Xtt) ] admissible for (p (x),x{t)) and with 
direction 

\x,t) = (m(x,t),U(x,t),0) 

such that 

\m(x, t)\ > f (po(x) X (t) - \m(x, t)\ 2 ) . (6.2) 

Since z := (m,U) and K^° x are uniformly continuous in (x,t), there 
exists an e > such that for any (x,t),(xo,to) G with |x — xo| + 
|t — to| < £ ) we have 

(z(x,t),g (^)) ± {m(x ,t ),U(xo,to),0) C hinti^. (6.3) 

Step 2. Fix (xo,to) £ & f° r the moment. Now, let < (j) ro < 1 
be a smooth cutoff function on with support contained in a ball 
B ro (xo,to) C ^ for some ro > 0, identically 1 on B ro / 2 (xo, to) an d 
strictly less than 1 outside. Thanks to Proposition 15.11 and the iden- 
tification (m,U,q) — > M, for the admissible line segment o~r XOt t )i there 
exist an operator Aq and a direction r/Q G R n+1 not parallel to e n +i, such 
that for any k G N 

/ cos(fcr/ • . 

V P J = (x °'* o) sm (^° ' 

and such that the pair (fhk,o, Uk,o) defined by 

(fn kfi ,U kfi )(x,t) := A [4> ro {x, t)k~ 3 cos(kr] • (x,t))] 

satisfies (|3.ip with q = 0. Note that (m^o, Utfl) is supported in the ball 
B ro {xQ,to) and that 

(m fcj0 , C/ fc)0 ) - 4> ro (rn(x , t ), U (x , t )) sin(kr]o ■ (x, t)) 

< const (A ,rj , ||^»0 Ilea) T ( 6 - 4 ) 

since Aq is a linear differential operator of homogeneous degree 3. Fur- 
thermore, for all (x,t) G B rQ /2(xo,to), we have 

\m k fi(x,t)\ 2 = |m(x ,to)| 2 sin 2 (A:r/o ■ (x,t)). 

Since rjQ G M n+1 is not parallel to e n +i, from Lemma [5.21 we can see that 

1 2 1 1 / 1 — / , m2 



lim / \m kfi (x,t)\Ux = - \m(x ,t Q )\ dx 

*"►«> J B ro/2 (x ,t ) 1 J B ro/2 (x ,t ) 
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uniformly in t. In particular, using (|6.2p . we obtain 



lim / \fhk q(x, t)\ 2 dxdt > 

JB ro/2 (x ,t ) 

F 2 ( 2 \2 

7^ — 7 — \ /. x [Po(x )x(t )) - \m(x ,t )\ I \B r/2 (x ,t )\ . (6.5) 



5iep 5. Next, observe that since m is uniformly continuous, there 
exists an r > such that for any r < f there exists a finite family of 
pairwise disjoint balls B rj (xj,tj) C with j-j < f such that 

2 \2 

/°o(^)x(^) ~~ \m(x, t)\ I dxdt < 



2 ^2(po(xj)x(tj) -\m(xj,tj)f) \B r .{xj,tj)\. (6.6) 



Fix s > with s < min{f, e} and choose a finite family of pairwise 
disjoint balls B rj (xj,tj) C with radii rj < s such that (|6.6p holds. In 
each ball B2 rj (xj,tj) we apply the construction of S'iep 5 to obtain, for 
every k G N, a pair (m fc j, t/jtj). 

Final step. Letting (m^, [/&) be as in (|6.ip . we observe that the sum 
therein consists of finitely many terms. Therefore from (|6.3p and (|6.4p 
we deduce that there exists ko G N such that 

G Xq for all k > fco- (6-7) 

Moreover, owing to ()6.5|) and (|6.6p we can write 

lim / |mfc(x, t) — m(x, t)\ 2 dxdt = lim / \fhk j(x, t)\ 2 dxdt 



j 

2 



>C (p (x)x(t) - \m(x,t)\ 2 ^ dxdt. (61 



Since rrik -4 m, due to (|6.8p we have 
liminf H^/cll^m = WmW^ + liminf \\rrik — Tn\\ 2 

k— >oo ^ ' fc— >oo 

/°o(2 ; )x(^) ~~ \m(x, t)\ 2 } dxdt, (6.9) 
which gives fl05j) with /3 = /3(n) = P{F(n)). □ 

7. Construction of suitable initial data 

In this section we show the existence of a subsolution in the sense 
of Definition 14.21 Since the subsolution we aim to construct has to be 
space-periodic, it will be enough to work on the building brick Q and 
then extend the costruction periodically to W l . 
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The idea to work in the space-periodic setting has been recently adopted 
by Wiedemann |22| in order to construct global solutions to the incom- 
pressible Euler equations. 

Proposition 7.1. Let po £ Cp(Q;R + ) be a given density function as 
in Proposition \4-l\ and let T be any given positive time. Then, there 
exist a smooth function % : R — » R + , a continuous periodic matrix field 
U : R n -> S$ and a function q £ C^R; C£(K n )) suc/t </itrf 



div x [/ + V x (f = on R"xR (7.1) 



and 



e(p (£), 0, t/(x)) < — /or a// fx, t) £ R n x [0, T[ (7.2) 
n 

t) = p(po(x)) + — /or all x £ R n x R. (7.3) 
n 

Proof. [Proposition 17.1] Let us define [/ componentwise by its Fourier 
transform as follows: 

Uij(k) := ( nk t\ u X {p7k)) if i * j, 
\{n-l)\k\ J 

Uu(k) := ( ^~^L )pMb)- ( 7 - 4 ) 
\{n-l) \k\ J 

for every k ^ 0, and U(0) = 0. Clearly C/jj thus defined is symmetric and 
trace-free. Moreover, since p(po) £ Cp(R n ), standard elliptic regularity 
arguments allow us to conclude that U is a continuous periodic matrix 
field. Next, notice that 

e(p (x),0,U(x)) = X m ax(-U) =A (7.5) 

oo oo 

for some positive constant A. Therefore, we can choose any smooth 
function x on R such that x > n ^ on [0,7"] in order to ensure (|7.2p . 
Now, let q be defined exactly as in (|7.3|) for the choice of x just done. It 
remains to show that (|7.ip holds. In light of (|7.3p . we can write equation 
(|7,ip in Fourier space as 

n ^ 

fcjt/ij = feiP(po) (7.6) 

J'=l 



for A; £ Z n . It is easy to check that [/ as defined by (|7.4p solves (|7.6p 
and hence U and q satisfy (|T. 1[) 

□ 

Remark 7.2. We noie f/iaf the Holder continuity of pq would be enough 
to argue as in the previous proof in order to infer the continuity ofU. 
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Proposition 7.3. Let p$ G Cp(Q;R + ) be a given density function as in 
Proposition \4-l\ and let T be any given positive time. There exist triples 
(m, U, q) solving (13. ip distributionally on W 1 x M enjoying the following 
properties: 

(m,U,q) is continuous in M n x (M\{0}) and me C(R; H w (R n )), 

(7.7) 

U(-,t) = U(-) fort = -T,T (7.8) 
and 

q( x ) = p(p (x)) + ^ for all (x, t) G R n x R, (7.9) 

e(p (x),m(x,t),U(x,t)) < ^ /or all (x,t) El"x ([-T, 0[n]0, T]). 

n 

(7.10) 

Moreover 

|m(x,0)| 2 = po(x) X (0) a.e. mR™. (7.11) 

Proof. [Proposition 17. 3| We first choose q := q given by Proposition 17. 11 
This choice already yields (I7.9p . 

Now, in analogy with Definition 14.21 we consider the space Xq defined 
as the set of continuous vector fields m : M n x] — T, T[— > W l in C°(] — 
T,T[;Cp(Q)) to which there exists a continuous space-periodic matrix 
field U : R n x] — T, T[— > Sft such that 

div x m = 0, 

d t m + div x U + V x q = 0, (7.12) 



supp(m) cQx [-T/2,T/2[ (7.13) 
[/■(., t) = U(.) for all t G [-T,T[\[-T/2,T/2] (7.14) 



and 



e(p (a;),m(x,t),Z7(a;,t)) < ^ for all (x,t) G R n x]-T,T[. (7.15) 

n 

As in Section 4.1, X$ consists of functions m :] — T, T[—> H taking values 
in a bounded set B C H . On B the weak topology of L 2 is metrizable, 
and correspondingly we find a metric d on C(] — T,T[;B) inducing the 
topology of C(]-T,T[;H W (Q)). 

Next we note that with minor modifications the proof of Lemma 14.51 
leads to the following claim: 
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Claim: Let Qq C Q be given. Let m G Xq with associated matrix 
field U and let a > such that 



|m(x,0)| 2 -(p (x)x(0)) 



da; < —a 



Then, for any 5 > there exists a sequence m k E Xq with associated 
smooth matrix field U k such that 

supp(m fc - m,U k -U) C Qq x [-5, <5], 

d 



and 



liminf / |mfc(x, 0)| 2 > / \m(x,0)\ 2 dx + f3a 



?o ■'Qo 
Fix an exhausting sequence of bounded open subsets Q k C (Jfc+i C Q, 
each compactly contained in Q, and such that |Qfc+i\Q/c| < 2 _fc . Let also 
7 e be a standard mollifying kernel in M n (the unusual notation 7 e for the 
standard mollifying kernel is aimed at avoiding confusion between it and 
the density function). Using the claim above we construct inductively a 
sequence of momentum vector fields m k G Xq, associated matrix fields 
U k and a sequence of numbers r] k < 2~ k as follows. 

First of all let mi = 0. Ui(x, t) = U(x) for all (x, t) G ~R n+1 and having 
obtained {m\,U\), (m k ,U k ), we choose % < 2~ k in such 

a way that 

\\m k -m k * j m \\ L2 < 2" k . (7.16) 

Then, we set 

«fc = - / [\m k (x, 0)| 2 - po(x)x(0))]dx. 
JQk 

Note that (17. 15ft ensures a k > 0. Then, we apply the claim with Q k , 
a = a k and 6 = 2~ k T to obtain m k+ i G Xo an d associated smooth 
matrix field Ut+i such that 



supp(m fc+1 - m k , U k+1 - U k ) C Q k x [-2~ fe T, 2- fe T], (7.17) 
d(m k+1 ,m k ) < 2~ k , (7.18) 

\m k+1 (x,0)\ 2 dx > / |m fc (x, 0)| 2 dx + /3o|. (7.19) 

Since d induces the topology of C(] — T, T[; H W (Q)) we can also require 
that 

\\{m k - m k+1 ) * 7 % || i2(n) < 2~ fc for all j < k for t = 0. (7.20) 
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From ([7]) we infer the existence of a function m E C(] — T, T[, H W (Q)) 
such that 

d 

m k — > m. 

Besides, ([7.17P implies that for any compact subset S of Qx] — T, 0[U]0, T{ 
there exists ko such that (mk, Uk)\s = (^fe j ^fco)ls f° r au k > ko. Hence 
(m k , U k ) converges in C® oc (Qx] — T, 0[U]0, T[) to a continuous pair (m, C7) 
solving equations (1712]) in M n x] - T, 0[U]0, T[ and such that fFT^ - fTT^ 
hold. In order to conclude, we show that also ()7.1ip holds for fn. 

As first, we observe that (|7.19|) yields 

ttfe+i < a k - f3a\ + \Qk+i\Qk\ < "fe ~ /3a 2 k + 2~ fc , 
from which we deduce that 

ak — > as k — > oo. 
This, together with the following inequality 



> 



\m k (x, 



po(x) X (0) dx > -(a k +C\Q\Q k \) > -(a k +CT 



implies that 



lim 



|m fc (a;,0)| 2 - po(x)x(0) 



dx = 0. 



(7.21) 



On the other hand, owing to (|7.16p and (|7.20p . we can write for t = 
and for every k 

\\m k — m\\ L 2 

< \\m k - m k * 7 % || L2 + \\m k * j m - m * i m \\ L2 + \\m * -y Vk - m\\ L2 



< 2~ k + \\m k +j * l m - rrik+j+i * l Vk \\ L 2 + 2" 

3=0 

< 2~ (fc_2) . 



(7.22) 



Finally, (|7.22p implies that m k (-,0) — > m(-,0) strongly in H(Q) as 
k — >■ oo, which together with (17.21 p gives 

2 



\m(x, I 



Po( x )x(ty f° r almost every x € 



□ 



8. Proof of the main Theorems 

Proof. [Proof of Theorem 12. 1| Let T be any finite positive time and 
Po £ Cp(Q) be a given density function. Let also (m,U,q) be as in 
Proposition 17.31 Then, define x(t) '■= x(*)> Qo( x ) '■= q( x )> 

mix, t) for t G [0, T] s 
™o(M) = _ 7 / ' J r7797 . (8-23) 
m(x, t — 21 ) tor t G [i , 2i J , 
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^t) = {^' t)f ° r * G[0 ' r] r , (8.24) 
V 7 | U(x,t-2T) forte [T,2T]. V ' 

For this choices, the quadruple (mo, Uo,qo, x) satisfies the assumptions 
of Proposition 14.11 Therefore, there exist infinitely many solutions m G 
C([0, 2T],H W (Q)) of ([23]) in M n x [0, 2T[ with density p , such that 

m(x, 0) = m(x, 0) = mix, 2T) for a.e. 

and 

| m (.,t)| 2 = p (-)x(0) for almost every (x,t) G M n x]0,2T[. (8.25) 

Since |mo(-,0)| 2 = po(")x(0) a - e - m as well, it is enough to define 
m°(x) = mo(x, 0) to satisfy also (|2.10p and hence conclude the proof. □ 

Proof. [Proof of Theorem \2.2\ Under the assumptions of Theorem \2.1\ 
we have proven the existence of a bounded initial momentum mP allowing 
for infinitely many solutions m G C([0, T];H W (Q)) of QT® on R™ x [0, T[ 
with density po- Moreover, the proof (see Proof of Proposition IT. 1 P 
showed that for any smooth function x '■ ^ — * ^ + with x > n ^ > the 
following holds 

\m(x,t)\ 2 = po(x)x(t) a.e. in R n x [0,T[, (8.26) 
\m°(x)\ 2 = p (x)x(0) a.e. inl" (8.27) 

Now, we claim that there exist constants C± , Ci > such that choos- 
ing the function x{t) > n ^ on [0; T[ among solutions of the following 
differential inequality 

X'(t) < -CiX 1/2 (t) - C 2X 3/ \t), (8.28) 

then the weak solutions (po,m) of (|2 . 6[) obtained in Theorem 12.11 will also 
satisfy the admissibility condition (|2.5p on W l x [0, T[. Of course, there 
is an issue of compatibility between the differential inequality (|8.28p and 
the condition x > this motivates the existence of a time T > 
defining the maximal time-interval in which the admissibility condition 
indeed holds. 

Let T be any finite positive time. As first, we aim to prove the claim. 
Since m G C([0, T]; H W (Q)) is divergence-free and fulfills (fB~26> (ET27I ) 
and po is time-independent, (|2.5|) reduces to the following inequality 

V(*> + m • V + + ^™ ■ V (-U < 0, (8.29) 

2 y Po(x) J 2 \Po{x)J 

intended in the sense of (space-periodic) distributions on W 1 x [0, T\. As 
Po G Cp(Q), there exists a constant Cq with po < Cq on M n , whence (see 

dS5SD - (iEZZD ) 



\m(x,t)\ < c Vx(t) a.e. on R n x [0,T[. (8.30) 
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Similarly we can find constants ci, C2 > with 




< ci a.e. in M ; 



(8.31) 



(8.32) 



As a conseguence of fj8.30j) - (18.32j) . (18. 29j) holds as soon as \ satisfies 



Therefore, by choosing C\ := 2ciCo and C2 '■= c 2 cq we can conclude the 
proof of the claim. 

Now, it remains to show the existence of a function x as m the claim, 
i.e. that both the differential inequality (I8.28P and the condition \ > n ^ 
can hold true on some suitable time-interval. To this aim, we can consider 
the equality in (|8.28p . couple it with the initial condition x(0) = xo f° r 
some constant xo > n ^ an d then solve the resulting Cauchy problem. 
For the obtained solution x> there exists a positive time T such that 
X{t) > n\ on [0,T[. 

Finally, applying the claim on the time-interval [0, T[ we conclude that 
the admissibility condition holds on K n x [0, T[ as desired. □ 

Proof. [Proof of Theorem 11.1] The proof of Theorem 11.11 strongly relies 
on Theorems l2.HI2.2i Given a continuously differentiable initial density 
p° we apply Theorems 12 . ltf2~2l for po(x) := p°(x) thus obtaining a positive 
time T and a bounded initial momentum m° allowing for infinitely many 
solutions m G C([0, T]; H W (Q)) of ([22]) on IR" x [0,T[ with density p° 
and such that the following holds 



for a suitable smooth function x '■ [0,T] — > Now, define p(x,t) = 

po (#) 1 [0 T[ • ^ ms snows that (|2.4p holds. To prove (|2,3p observe that 
p is independent of t and m is weakly divergence-free for almost every 
< t < T. Therefore, the pair (p, m) is a weak solution of (jl.ip with 
initial data (p°,m°). Finally, we can also prove (12 . 5|) : each solution 
obtained is also admissible. Indeed, for p(x,t) = po(x)lr yr(t), (|2.5p is 
ensured by Theorem 12.21 □ 



X'(t) < -2cic x 



/ 2 (t)-c 2 c oX 3/2 (t) on [0,T[. 



m(x,t)\ 2 = po(x)x(t) a.e. in R n x [0,T[, 
m°(x)| 2 = p (x)x(0) a.e. in M n , 



(8.33) 



(8.34) 
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